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ABSTRACT 

The edge domain decomposition method (edge-DDM) is developed as a very powerful extension of the domain 
decomposition methods presented in [1] for the rapid and accurate simulation of light scattering from advanced photomasks. 
The range of validity of the method is systematically evaluated and it is found to be accurate for wavelength-sized mask fea- 
tures with large vertical topography. The error associated primarily with neglecting comer effects is seen to concentrate at the 
extremities of the spectrum of propagating plane waves and is therefore filtered-out by the NA and the larger than one reduc- 
tion factors in typical projection printing systems. A possible algorithmic implementation of the method, that involves rigor- 
ous pre-calculation of the edge-diffraction of all types of edges present in the layout, is presented for the simulation of large 
and arbitrary 2D layouts of alt. PSMs. A speed-up factor of 172,800 (Isec vs. 2days) is shown in the simulation of simple 3um 
by 3um (4X) 2D layouts (isolated hole and isolated island) with accuracy better than 99% compared to the rigorous 3D simu- 
lation of the mask diffraction. An example of the edge-DDM applied on a large 12um by 16um 2D layout of a 0°/90°/270° alt. 
PSM, where the near fields are estimated in less than lmin, is also presented. Spectral matching of the edge-diffraction with 
more compact piecewise constant models is shown to lead to additional speed-ups of the edge-DDM. 

Keywords: domain decomposition methodVedgeidiffraction, phase-shift mask modeling, diffraction order 

1. INTRODUCTION 

Researchers have in the past investigated the decomposition of mask layouts into individual openings, and even into 
edges, for ideal thin masks under coherent illumination [6], but such techniques have never been applied to the rigorous solu- 
tion for the scattered field. Domain decomposition techniques, where a large electromagnetic problem is broken up in smaller 
pieces and the final solution is arrived at by synthesizing (field-stitching) the elemental solutions, have been recently proposed 
for the study of one-dimensional binary (phase only) diffractive optical elements [5], [8J. Another research group working on 
the same problem demonstrated how to create and use a perturbation model for binary edge-transitions based on the product of 
the ideal, sharp transition and the continuous field variations in the vicinity of the edge [4]. Independently, a similar technique 
to field-stitching has been invented for the simulation of two-dimensional layouts of advanced photomasks (alternating PSM, 
masks with OPC) that properly models interactions from neighboring apertures and furthermore takes advantage of the spec- 
tral properties of the diffracted fields to come up with a compact model for the edge-transitions [1], 

In this paper, a domain decomposition method based on edges is described, where the pre-calculated electromagnetic 
field from the diffraction of isolated edges is recycled in the synthesis of the near diffracted field from arbitrary two-dimen- 
sional diffracting geometries. Although at first glance it might be tempting to discard edge-decomposition techniques as inac- 
curate for masks with large vertical topography and features with small lateral dimensions, it will be shown that for a large set 
of practical situations they are very accurate. A key difference in this method from the one in [4] is that the sum of the complex 
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Figure 1. Decomposition of an opening into two edges 
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referred to as the edge domain decomposition method or edge-DDM. The application of the edge-DDM to a single isolated 
space and an isolated line are shown next, followed by a systematic evaluation of the limits of applicability of the method. 

2.1. Decomposition of a space into two edges 

Consider the isolated mask space (opening) shown in Figure 2(a). The amplitude of the complex Ey-field under Ey 
(TE) illumination is shown in (c) and the scattered field across the observation plane is shown in (e). According to Figure 1, 
the field across the observation plane can be alternatively obtained by using the scattered field from an edge that has the same 
vertical profile as the space, as shown in Figure 2(b) and (d). The resulting field from this edge-DDM is overlaid on the plot of 
Figure 2(e) along with the field from the edge scattering. The normalized mean square error 1 (NMSE) between the r-mask and 
the edge-DDM is 0. 19%. In Figure 2(f) the spectra of the r-mask and edge-DDM are compared and it is seen that they differ by 
0.13% in a NMSE sense if the whole spectrum of propagating waves is considered and by merely 0.007% within the collection 
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Figure 2. Edge-DDM applied on an isolated space 

Geometry of mask with isolated space (a) and edge (b). Note that the vertical profile of the edge is the same as that of 
the isolated space. Amplitude of complex field from rigorous TEMPEST simulation for space (c) and edge (d). The 
true scattered field (r-mask) across the observation plane in (c) is compared with the edge-DDM (synthesized field) 
in (e) and their spectra are compared in (f). The normalized mean square errors (NMSE) are indicated on the plots. 
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1. The normalized mean square error is defined by: NMSE = J|E r DE s | dS/ j]E r | dS , where E, is the "true" field resulting from 
a rigorous simulation of the complete mask geometry and E s is the synthesized field via the edge-DDM. 
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ability of a projection system witb NA=0.8 and R=4. Since only the portion of the spectrum that is collected by the projection 
system contributes to the image formation, the exceptional accuracy of the edge-DDM in this regime is more than adequate for 
accurate image simulations. 

2.2. Decomposition of a line into two edges 

Similar steps can be followed for the isolated line shown in Figure 3(a) uhumnated with E, (TM) polarized light 
Note that it is not typical to have a line as shown here, where in both sides the glass has been etched. Normally only one side 
of the line is phase-shifted, but this example simply aims to introduce the principle of the edge-DDM. A line with one side 
phase-slnfted and the other left intact can still be decomposed via the edge-DDM using the diftraction fields from the two dif- 
ferent edges. The amplitude of the complex E, field everywhere in the domain of Figure 3(a) is shown in (c) and for the edge 
of Figure 3(b) in (d) respectively. The true (r-mask) scattered field is compared with the synthesized field with edge-DDM in 
(e) and the spectra are compared in (f). The error (NMSE) is seen to be very small in both the near field and spectra plots and 
specifically the error of the through-the-lens (TTL) spectrum is just 1.5E-5. 
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Figure 3. Edge-DDM applied on an isolated line 



Geometry ^of mask .with isolated line (a) and edge (b). Note that the vertical profile of the edge is the same as that of 
^ ST ^ Phh ! dC ° f COmpleX fidd fr0m ^rous TEMPEST simulation for line (c) and edge (d) The 

X5S^ S2ET* aCr0SS ^ e ^T 10n P,ane in (C) is ^ *e edge-DDM (synthesized Afield) 
in (e) and their spectra are compared in (f). The normalized mean square errors (NMSE) are indicated on the plots 



c 



c 



23. Limits of the edge-DDM 

In the examples of Figure 2 and Figure 3 the edge-DDM is seen to be very accurate compared with the simulation of 

the exact mask structure. Although the vertical mask topography ( 1 80° * 170nm « 1 X ) is comparable to the wavelength, the 
fact that the lateral dimensions of the space and the line are relatively large compared with the wavelength 
(CD = 400nm > 2 x 193nm = 2 x X ) renders the edge-DDM accurate. One expects that as the lateral sizes become smaller 
and the vertical mask topography even larger the edge-DDM will eventually break down. The limits of the edge-DDM were 
investigated through exhaustive simulation of isolated spaces and lines on the mask with various etched depths and amounts of 
undercut. The sidewall angles were always kept vertical and no rounding at the bottom of the etched wells was included. The 
mask CD was varied in the set {2\im, l*im, 0.6nm, 0.4nm, 0.2^m(-lA;), O.ljim}, the etched depth in {0°, 90°, 180°, 270°, 
360°} and the amount of underetch in {Onm, 25nm, 50nm, 75nm, lOOnm). All 2 (spaceAine) by 2 (TE/TM illumination) by 6 
(CD sizes) by 5 (etched depths) by 5 (underetch) = 600 cases were run and the rigorous (r-mask) scattered field was automati- 
cally compared with that obtained from edge-DDM. The results are spectacular! The edge-DDM can readily achieve better 
than 1% accuracy (in a NMSE sense) in the near field for CD=400nm (4X) or larger regardless of depth, underetch and polar- 
ization. The accuracy in the TTL-spectrum is sufficient in those cases even for inspection (R=l) simulations. For lines, the 
edge-DDM maintains excellent accuracy in the TTL-spectrum with R=4 down to a mask CD=100nm (~0.5X) with 90° of 
etched depth and Onm of undercut, or mask CD=200nm (~1 A.) with 270° of etched depth and 25 nm of undercut For spaces, the 
accuracy of the method is acceptable at least down to mask CD=200nm (~1X) with 270° of etched depth and 50nm of 
underetch. 

Current state-of-the-art photomasks and processes utilize sub-resolution assist features that are never less than one 
wavelength in size on the mask. This clearly implies that the edge-DDM has the potential to accurately simulate the most 
advanced photomask technologies (OPC, alt. PSM, OPC on alt. PSM). As with the DDM of [1], the true power of the method 
lies on its application in 2D mask layouts. However, one key difference with the DDM of [1] is the following: Only one rigor- 
ous 2D simulation of an edge suffices for the scattered field reconstruction of any arbitrary size on the mask! The DDM of [1] 
requires a separate simulation of every different mask size. Next, the application of the edge-DDM in 2D layouts is developed. 

3. THE EDGE-DDM APPLIED IN 2D LAYOUTS 

The application of the edge-DDM in 2D layouts is straightforward. A mask layout comprises of a large number of 
edges positioned at different locations and having different orientations. No matter how complicated the mask technology is, 
there is usually only a small number of different types of edges present in the layout. For example, a single exposure 0°/90°/ 
270° alt. PSM has five types of edges: i) Cr-layer/0° edge, ii) Cr-layer/90° edge, iii) Cr-layer/270° edge, iv) 0°/90° edge and 
0°/270° edge. Depending on the orientation of each edge in the layout it "sees" and responds to the incident field differently. If 
the illuminating field is a TE (Ey) normally incident plane wave and the edge is oriented along the y-axis it "sees" TE illumi- 
nation, but if it is oriented along the x-axis it "sees" TM illumination. This is illustrated in Figure 4. The response of each edge 
to its respective illumination is taken into account rigorously from a pre-stored 2D edge diffraction simulation. However, cor- 
ner effects are effectively ignored, since the finite extent of each edge is not rigorously taken into account. Instead, the scat- 
tered field at the end points of every edge is abruptly terminated (truncated) in a perfect square- wave fashion to the field value 
of the k-mask model. Comer effects in typical imaging situations are mapped at the extremities of the spectra and do not con- 
tribute to the image formation. Due to this observation, which will be shortly proven valid in examples to follow, the edge- 



c 



c 



DDM applied in 2D layouts is very successful in accurately capturing the true electromagnetic behavior of edges in a very 
rapid manner, while the insignificant comer effects are safely ignored. 
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Figure 4. Edge-DDM applied on an arbitrary Manhattan-type polygon 
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Slightly different embodiments of the edge-DDM applied in 2D layouts can be devised. One that is very powerful 

m alt. PSMs are presented for completeness. Next, application of the edge-DDM in simple 2D mask layouts is presented and 
the results of the method are compared with fully rigorous r-mask models in order to establish the validity of this approach 
Finally a large portion of a real layout from a single exposure 0°/90°/270° alt. PSM that is too big for 3D rigorous simulation 
is simulated via the edge-DDM and the aerial image is compared with the simple k-mask model. 

3.1. Algorithmic implementation of the edge-DDM 

First, every distinct edge that is present on the layout is pre-simulated with all possible illumination directions (field 
pobnzauons) that are required, based on the orientations of the edge encountered within the layout. For simplicity only Man- 
hattan-Type layouts are considered, where the edges are oriented along either the x- or the y-axis and require electromlgnetic 
simulation under both parallel and perpendicular to the ^,--^1^.^,..^!^,^^ 
much more computationally .aborious is straightforward. Subsequently, each distinct mask layer is broken up into a set of 
mutuaUydisjomt^on-overlapping) rectangle, For t^^l^^^^^^^*^ 
layer (typically Cr-based) and the clear (uncovered) layer, whereas a more advanced 0°/90°/270° three-phase alt PSM consists 
of four layers, namely the absorption layer and the 0° clear, 90<> clear and 270° clear layers. Next, looping through all rectan- 
gles o each layer, the y of each one of the four edges of the rectangle is determmed based on what its neighbo^g layer is 
Trivial edges (when the neighboring rectangle is of the same type) are discarded. It is also possible that only part of a Lgle 
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edge from a rectangle is neighboring with one layer and the rest is neighboring with a different layer (or different layers). In 
that case the edge is broken up in multiple edges, in a way that tracks the neighboring layer. What is left at the end of this pro- 
cess is a set of all non-trivial edges that are present in the layout (location, size and orientation) and their type, te. which are 
the two mask layers on each side of the edge. Given the mask layout, the k-mask (Kirchhoff-mask) model is readily available. 
Finally, the respective (complex) difference of the true edge scattering from the ideal and shaip k-mask model is added to all 
non-trivial edges accordingly. This revised mask model that results from the edge-DDM includes accurate information about 
the electromagnetic scattering from the edges and it can be thought of as a quasi-rigorous mask model (qr-mask) in the sense 
that only 2D rigorous mask simulations are used to approximate a 3D electromagnetic problem. 

The above algorithm was implemented in the MATLAB environment. The part that deals with the polygon and edge 
extraction can probably be judged as rudimentary and most likely cannot compete in speed with more sophisticated implemen- 
tations in production caliber CAD software, but it suffices for the purposes of this work, that is, to demonstrate the proof of 
concept of the edge-DDM and the feasibility of including very accurate edge scattering information in rapid aerial image sim- 
ulations. 

3.2. Dependence of edge scattering on profile and polarization 

Examining the departure of the scattered field from an edge from its ideal, sharp edge transition and how it depends 
on the profile of the edge and the incident field polarization can provide great intuition of the physical mec h a nism s involved. 

Figure 5 shows the amplitude of the (near) scattered field from an edge with various profiles. The ideal, sharp edge 
transition is also shown on every graph for comparison. In Figure 5(a) and (b) the edge scattering for different etched depths 
(0°, 180° and 360°) and vertical sidewalk with no undercut is shown for TE (Ey) and TM (EJ illumination respectively. The 
effect of undercut is depicted in Figure 5(c) and (d), where the edge scattering for different amounts of undercut (Onm, 25nm, 
50nm and 75nm) and vertical sidewalk with 180° etched depth is shown for TE and TM illumination respectively. It is inter- 
esting to note that increasing the etched depth is "making the edge response slower", meaning that the near scattered field from 
the edge requires more distance to reach its "on" value from its "off* value, although thk effect is more pronounced for the 
first increment from 0° to 180° rather than from 180° to 360°. Similarly, increasing the amount of undercut appears to "make 
the edge response faster" and furthermore this behavior appears to be linear, since for every increment of the undercut the edge 
transition gains a constant distance in its "off-on" behavior. 
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Figure 5. Dependence of edge scattering on profile and polarization 
Normalized amplitude of the near field across the observation plane for an edge with 0° 180° and 360° of e trh*H 
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It is instructive to isolate the difference of the true edge scattering from the ideal, sharp edge transition. This is done 
in Figure 6, which depicts the amplitude of this (complex) difference for all respective plots of Figure 5. Observe that, indeed, 
by increasing the etched depth the difference has a larger spatial extent, whereas increasing the amount of undercut confines 
the spatial extent of the difference. This is crucial, because it is exactly the spatial extent of the difference that affects its spec- 
tral distribution. If the difference is well confined, its spectral distribution is spread out to higher spatial frequencies and the 
lower frequency content is small. In such a case the complex difference can be neglected and the approximation of the ideal, 
sharp transition is adequate. However, when the difference has a larger spatial extent then its lower frequency content becomes 
significant and it cannot be neglected. This is shown in Figure 7, which depicts the magnitude of the spectra of all respective 
difference distributions of Figure 6. Observe that the through-the-lens (TTL) spectrum of the difference attains larger values 
when the etched depth increases and smaller values when the amount of undercut increases, as expected from Figure 5 and 
Figure 6. The normalization used in spectra plots is that a constant (DC) value of 1 (clear field) would produce a Dirac 8-func- 

tion at zero frequency with amplitude of 19.4 = Jj\Lpz 0 * WACl). 
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Figure 6. Complex difference of true edge scattering from ideal, sharp transition 

The amplitude of the complex difference of the true electromagnetic field scattered from an edge from the ideal, sharp 
transition (step-function) is shown in (a)-(d) for all edge profiles and illumination polarizations of Figure 5 respec- 
tively, (a) and (c) are for TE polarization and (b) and (d) are for TM polarization. The spatial extent of the difference 
grows with increasing etched depths, as seen in (a) and (b), and shrinks with increasing undercuts, as seen in (d) for 
TM polarization. The behavior in (c) appears somewhat erratic, but if the phase is also taken into account, it is shown 
in Figure 7 that this is also the case for TE polarization. 
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Figure 7. Spectrum of complex difference 

The spectra of the complex differences shown in Figure 6 are depicted here respectively, (a) and (c) are for TE polar- 
ization and (b) and (d) are for TM polarization. Observe that the through-the-lens (TTL) parts of these spectra, which 
l ! thog ^ p 1 hlc ""S^* situations NA-0.7-0.8) extent up to ~ +/- 0.008 (nm)" 1 are smaller for smaller 
etched depths and larger undercuts, regardless of polarization, as expected from the observations on the plots of 
Figure 5 and Figure 6. r 

3.3. Application of edge-DDM in simple 2D layouts 

The application of the edge-DDM on an isolated square hole and an isolated square island (post) and comparison of 
the results with fully rigorous 3D mask simulations (r-masks) are presented here. The edge profile has a glass etched depth of 
180° and an undercut of 50nm. The dimensions of both the hole and the island are 400nm x 400nm (4X), corresponding to 
approximately 2x2 wavelengths, for X=193nm. 

Figure 8(a) and (b) depict the amplitude of the scattered field across the observation plane below the hole obtained 
via the edge-DDM and via fully rigorous 3D simulation, respectively. The amplitude of their difference (error of the edge- 
DDM) is shown in Figure 8(c). Locally the amplitude of the error is seen to reach a discouraging level of almost 50% of the 
clear field value. The normalized mean square error is however more contained and is only about 4.9% if calculated within a 
3pm by 3pm area or 3.6% if calculated within the 0.4pm by 0.4pm open area of the hole, as shown in Figure 8. However, nei- 
ther of these error levels is immediately relevant in the simulation of the image formation. What is rather more relevant is the 
error incurred on the dif&action orders that are collected by the imaging optics. If the spatial frequencies of the error depicted 
in Figure 8(c) are sufficiently high, this error will be discarded by the projection lens. This concept is shown in Figure 9, where 
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Figure 8. Edge-DDM applied on an isolated square hole (scattered field) 

(a) Synthesized with the edge-DDM and (b) rigorously calculated scattered field (amplitude), across the observation 
plane below a 400nm x 400nm isolated hole. A 1 80° etch and 50nm undercut has been applied. The amplitude of the 
error (difference of r-mask from edge-DDM) is shown in (c). Locally the amplitude of the error reaches the 50% 
level, but the normalized mean square error does not exceed 5%. 

the spectra (magnitude) of the scattered fields from edge-DDM and from rigorous 3D simulation are shown in (a) and (b) 

respectively. The marginal circles indicating the propagating plane waves (with |k| <, 2n/X ) and the collected by the imaging 

system plane waves (with |k| £ 2tiNA( 1 + a)/XR ), for NA=0.75, R=4 and o=0.3, are indicated on these plots. The magni- 
tude of the error in the spectrum of the field obtained from edge-DDM is shown in Figure 9(c). Observe that most of the error 
is concentrated in the higher spatial frequencies for |k x | approaching 0.03 and only small levels of error exist in the collected 
orders. The normalized mean square errors of all propagated and collected orders are 2.2% and 2E-4 respectively. 
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Figure 9. Edge-DDM applied on an isolated square hole (spectrum) 

Magnitude of the plane wave spectra of the near field plots in Figure 8. Spectrum of the synthesized (edge-DDM) 
field in (a) and of the rigorous field (from 3D simulation) in (b). The spectrum of the error is depicted in (c). The error 
level within the collected orders circle is seen to be an insignificant 0.02%. 
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Figure 10. Edge-DDM applied on an isolated square island (scattered field) 

tude of the error (difference of r-mask fiWed^nni^ ;! a w fT !. derCUt 1,38 been "n* 84 1,16 ■«**- 
the 50% .eve,, but the no^L^^^sl^ «* - ^de of the error reaches 



(a) 



Propagated orders 



o 

2k o 



o I 



SI 



\ Collected orders / 



>.03 -0.01 0.01 
k,(nm)> 



0.03 




(C) 



Si 



O . 

o 



9 



\ Collected orders / 



0.13 



21 



-001 0.01 
Mnm)-' 



'-0.03 -0.01 0.01 0.03 
Mnm)-' 

NMSE(propagated orders) = 4% 
NMSE(collected orders) = 0.16% 
Figure 1 1 . Edge-DDM applied on an isolated square island (spectrum) 

level within the collect orders circieTonly EST* W ° f *" ^ is de P icted » «* ™* error 



7 



C 



It is worth mentioning that each fully rigorous mask simulation that was performed with TEMPEST used approxi- 
mately 1.8Gbytes of memory and required almost 2 days on a 550Mhz CPU. On the other hand, once all of the required 2D 
edge simulation results are available the scattered field can be obtained with the edge-DDM in less than lsec! It should be 
stressed once again here that the simulation of the isolated edge diffractions is an off-line procedure. It only needs to be per- 
formed once and the edge-diffraction results are recycled as many times as are necessary for the reconstruction of the scattered 
fields of arbitrary layouts. In the above examples, it took approximately 5min for each rigorous, isolated-edge 2D simulation. 

3.4. Example of edge-DDM on a large, arbitrary layout of a 0°/90°/270° alt PSM 

The edge-DDM is used here to generate the scattered field from a three level (0°/90 o /270°) alt. PSM. A portion of the 
arbitrary target layout is shown in Figure 12 before and after the phase assignment. The rectangle and edge extraction accord- 
ing to the algorithm presented in Section 3.1 is performed within the 3^mby 4^im (IX) coherence window shown and the scat- 
tered field obtained with the edge-DDM inside the 12^im by 16^m (4X) coherence window is also shown in Figure 12 for TE 
and TM illuminations at k=193nm. Note that a fully rigorous 3D mask simulation of a 12|im by 16jjm (4X) layout would 
require close to 40Gbytes of memory and would only be feasible on a multi-CPU architecture. The scattered field is obtained 
with the edge-DDM implementation in less than lmin. This time includes the polygon and edge extraction and the field syn- 
thesis from pre-calculated edge diffractions. Therefore, the results of the edge-DDM cannot be immediately compared with a 
fully rigorous solution, but the accuracy of the edge-DDM within the imaging system's collection ability can reveal the inaccu- 
racies of the ideal k-mask, where all transitions are assumed perfectly sharp. This is done in Figure 13, where the images (at 

Target layer After phase assignment Rectangle and edge 




Figure 12. Example: Edge-DDM on a large 2D mask layout (algorithm and near fields) 

The target mask layer and the mask after assignment of an arbitrary layout are shown. A three phase-level alt. PSM is 
assumed (0°/90°/270°). The rectangle and edge extraction within the 3\xm by 4|im (IX) coherence window for the 
various polygons is also shown. Finally, the scattered near field within the coherence window for both TE and TM 
(X=193nm) illumination is depicted. 
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Figure 13. Example: Edge-DDM on a laige 2D mask layout (images) 
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the 30% intensity level) obtained with the k-mask and the edge-DDM for two different undercuts (Onm and 50nm) are com- 
pared. The system parameters are X=193nm, NA=0.7, o=0.3 and R=4. The inaccuracies of the k-mask model are revealed at 
vanous locations in the image as seen in Figure 13. The effect on the image of changing the undercut from Onm to 50nm can 
be quickly and accurately evaluated. Observe the subtle differences evident in the intensity contours of the Onm and 50nm 
masks of Figure 13. 

4. SPEEDING-UP FURTHER THE EDGE-DDM 

Although the edge-DDM presented earlier achieves excellent accuracy and tremendous speed-ups as compared to a 
fully rigorous mask simulation, there is still one more parameter to exploit in order to further speed-up the method As was 
shown in [1], where the matched Kixchhoff mask model (mk-mask) was developed, there is no need to carry all the details of 
the true field scattering as far as the imaging process is concerned (with larger than 1 reduction factors), since the higher spa- 
tial frequencies are not used. If a lumped parameter model can be devised that closely matches the through-the-lens (TTL) 
spectrum of the scattered field then the task of an accurate field representation is accomplished. In [1] a rect-function was seen 
to very successfully mimic the TTL spectra of scattered fields from vanous size and profile openings. A similar principle can 
be apphed here with the purpose of matching the TTL spectra of edge-diffracted fields. Assuming a restriction to piecewise 
constant models, the scattered field from an isolated edge under TE or TM Ulumination can be approximated using the multi- 
step functions shown in Figure 14. Although not shown in Figure 14, the amplitude, phase and size of all steps are optimized 
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Figure 14. Matched bandwidth approximation to edge diffraction 

Instead of approximating the true edge diffraction by a single step function (geometrical or KirchhofF approximation) 
a multi-step piecewise constant function is employed. The amplitudes, phases and sizes of all steps are optimized to 
achieve the best spectral matching of the collected orders with the orders resulting from the continuous edge diffracted 
field. The matched bandwidth approximations to the field diffracted by a 1 80° edge with 50nm undercut under TE and 
TM illuminations are shown in (a) and (b) respectively. 

to achieve the best spectral matching of the collected (TTL) orders with the orders resulting from the continuous edge-dif- 
fracted field. Note that the key factor for a good spectral match is to capture the finite extent of the "off-on" edge transition. A 
piecewise constant function is not the best choice to achieve this task, but on the other hand a piecewise constant function 
offers a good lumped parameter approach and it lends itself nicely as an extension of the unmatched Kirchhoff approach, 
where all edge transitions are assumed perfectly sharp. 

If the matched bandwidth edge-DDM is applied on the 400nm by 400nm (4X) isolated hole of Figure 8 the scattered 
field across the observation plane can be represented in the compact way shown in Figure 15(a). Recall that the edges parallel 
to the y-axis "see" TE polarization and the edges parallel to the x-axis TM polarization. The multi-step matched BW approxi- 
mations to the edge diffraction shown in Figure 14 for TE and TM illuminations intermingle in the x- and y-directions to pro- 
duce the multi-color mask of Figure 15(a), where each color represents a different amplitude and phase region of the mask. 
The spectrum of the error of this approach is depicted in Figure 15(b). Note that compared to the error shown in Figure 9(c) 
this error has increased by almost four- fold. Nevertheless, the error in the TTL spectrum still remains at low enough levels for 
accurate imaging simulations. A better approximating function, other than the piecewise constant step function, that is more 
capable of capturing the true edge-diffraction profile (such as linear or higher order polynomial) should result in better spectral 
matching and consequendy more accurate imaging simulations. 
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Figure 15. The matched BW edge-DDM applied to the hole of Figure 8 

(a) Representation of the scattered field below a 400nm by 400nm hole (etched depth 1 80°, undercut=50nm) using the 
matched bandwidth edge-DDM, as shown in Figure 14. Each color represents a different amplitude and phase region 
of the mask. The spectrum of the error of the representation shown in (a) compared to the true solution of Figure 8(b) 
is depicted in (b). The error level within the collected orders circle is seen to increase from 0.02% to 0.086%. 



5. CONCLUSIONS 

This paper has dealt with the development of the edge-DDM (edge domain decomposition method). This method is 
an extension of the domain decomposition method proposed and developed in [1] that introduces tremendous versatility, since 
only a small number of isolated edge diffraction simulations is shown to contain all the necessary information for the synthesis 
of the scattered field from arbitrary 2D mask layouts and subsequent accurate imaging simulations. The limits of this method 
are reached when the mask features are smaller than a wavelength in size and the vertical mask topography is large. Through a 
systematic process it was detennined that features as small as at least 200nm bearing 270° of phase-wells can be accurately 
decomposed via the edge-DDM, at A=193nm. The method was tested for simple 2D layouts, where rigorous mask simulations 
are possible. Excellent accuracy accompanied by speed-up factors of 172,800 (lsec vs. 2days) were demonstrated. The accu- 
racy of the method was attributed to the fact that on one hand the edge diffraction phenomena are modeled rigorously through 
2D edge-diffraction simulations and that on the other hand the errors incurred during the synthesis (primarily near the comers) 
are mapped at the extremities of the spectrum of propagating waves and do not contribute to the image formation. It was 
emphasized that all necessary rigorous 2D simulations of the diffraction from isolated edge-profiles is performed off-line and 
recycled for the diffraction calculation of arbitrary layouts. If needed, the accuracy of the edge-profile diffraction simulations 
can be pushed to extreme limits, since it is performed only once. This can aid for example in the correct simulation of unusual 
edge-profiles whose geometrical details would require excessive discretization of the domain. It was also shown that the exact 
details of the edge-diffraction are not necessary for accurate imaging simulations (with larger than 1 reduction factors) and a 
piecewise constant, multi-step edge transition model that matches the through-the-lens (TTL) spectrum was introduced The 
range of validity of the edge-DDM is expected to be appropriate for rapid and accurate evaluation of aerial images whenever 
speed is critical, as for example in full-chip OPC software and data-to-database comparisons in the inspection of alt. PSMs. 
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method*', in which the rect-CMTF, which is now a function of both x, y, is determined independently for the x- and y- direc- 
tions of the alt. PSM. A total of only four 2D rigorous EM simulations are required (shown in Figure 14(a) as cut-lines 1-4) 
that take just under 8min to run and utilize negligible memory (-3Mb). Note that in this example the cross-talk can be 
neglected (the 50nm of underetch virtually eliminated cross-talk between adjacent features) and the mask decomposition fol- 
lows Figure 3. 

The agreement of the "accurate scalar CMTF method" with the full 3D rigorous EM field simulation is spectacular. 
The image intensity contours at 30% of the clear field virtually overlap and the same is true for the aerial images at the cut-line 
that runs through the center of the contact holes, as shown in Figure 14(b) and 14(c). For comparison, the aerial image con- 
tours and at the same cut-line of the scalar CMTF without adjustment is also shown on these plots. 

Three-dimensional simulation results from other dark-field alt. PSM layouts (dense line/space pattern, corners) were 
also used to benchmark the "accurate scalar CMTF method" with equally successful outcomes. The ~400X speed-up achieved 
in the example presented would have been approximately halved for a discretization of 25 cells per A,. The speed-up resulting 
from this method depends on the size of the 2D mask layout, the discretization chosen for the rigorous EM field simulations 
and the number of rectangles that constitute the layout of the dark- field alt. PSM. 

6. CONCLUSIONS 

Rigorous electromagnetic simulation of the role of the mask in the imaging process has so far been limited to 2D 
problems and relatively small periodic 3D problems. Although this capability is enough to tackle a variety of advanced photo- 
mask design issues, it is unquestionably very restrictive. For example, rigorous simulation of a lO^m by lO^im area on an alt. 
PSM that contains random (non-periodic) structures, at X=193nm, would require on the order of ~3.6Gb of memory (discreti- 
zation of the ~14,000X, 3 domain with 20 cells per A.) and would run for many hours on a powerful CPU. Although the scalar 
theory can handle large mask domains, it has been discarded for advanced photomask technology simulations as being inaccu- 
rate. In this article, we have demonstrated a way to improve the accuracy of simulations with the scalar theory via the "accu- 
rate scalar CMTF method". This method involves using scalar CMTF (complex mask transmission function) models of single 
openings that are tuned to provide excellent spectral matching with the rigorously calculated CMTF inside a band of interest. 
A simple scalar CMTF model, namely the rect-CMTF, was presented and the effect of its three parameters on the spectrum 
was thoroughly investigated. The "accurate scalar CMTF method" is capable of better than 99% agreement with the vector 
theory and huge computational speed-ups, on the order of 200X or better, can be realized. Model-based CAD tools, that so far 
rely on scalar theory for aerial image calculations, can be tuned to be very accurate. Conceivably, accuracy in par with a rigor- 
ous 3D simulation becomes feasible for large - even full mask - areas. 

The way to decompose an alt. PSM into many constituents that properly account for the cross-talk between phase- 
wells was demonstrated, as an integral part of the "accurate scalar CMTF method". Through this decomposition, a much 
deeper understanding of the physical phenomena that take place at the corners and at the edges of phase-wells of alt. PSMs 
arose. Specifically, the comers and edges at the bottom of the phase- wells were identified as the sources of cross-talk between 
adjacent phase- wells and it was found that the amount of light due to cross-talk becomes significant for alt. PSM technologies 
deeper than 90°/270°. 
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